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Abstract 

Electroweak precision measurements can provide indirect information about the possible scale 
of supersymmetry already at the present level of accuracy. We update the present-day sensitivities 
of precision data using mt = 172.7 it 2.9 GeV for the experimental value of the top-quark mass, 
within the constrained minimal supersymmetric extension of the Standard Model (CMSSM), in 
which there are three independent soft supersymmetry-breaking parameters mi/2,mQ and ^o- In 
addition to Mw and sin^ ^eff > the analysis is based on {g — 2)^, BR(6 — > 57) and the lightest MSSM 
Higgs boson mass, M/j. Assuming initially that the lightest supersymmetric particle (LSP) is a 
neutralino, we display the CMSSM results as functions of mi/2) fixing mg so as to obtain the cold 
dark matter density allowed by WMAP and other cosmological data for specific values of tan /? 
and /i > 0. For a sample value of tan/3 we analyze how the global function would change 
following a possible future evolution of the experimental central value of mt and its error. In a 
second step, we extend the analysis to other constrained versions of the MSSM: the NUHM in 
which the soft supersymmetry-breaking contributions to the Higgs masses are independent and the 
Higgs mixing parameter fi and the pseudoscalar Higgs mass Ma become additional free parameters 
compared to the CMSSM, a VCMSSM in which the bilinear soft supersymmetry breaking parameter 
Bq = Aq — mo, and the GDM in which the LSP is the gravitino. In all scenarios we find indications 
for relatively light soft supersymmetry-breaking masses, offering good prospects for the LHC and 
the ILC, and in some cases also for the Tevatron. 
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1 Introduction 



We have recently analyzed the indications provided by current experimental data concerning 
the possible scale of supersymmetry [1,2] within the framework of the minimal supersymmet- 
ric extension of the Standard Model (MSSM) [3,4], assuming that the soft supersymmetry- 
brcaking scalar masses rriQ, gaugino masses mi/2 and tri-linear parameters Aq were each con- 
strained to be universal at the input GUT scale, with the gravitino heavy and the lightest 
supersymmetric particle (LSP) being the lightest neutralino Xi^ ^ framework often referred 
to as the constrained MSSM (CMSSM). However, this is not the only possible scenario for su- 
persymmetric phenomenology. For example, the soft supersymmetry-breaking scalar masses 
mo might not be universal, in particular those of the MSSM Higgs bosons, a framework 
we term the NUHM [5,6]. Alternatively, one may postulate supplementary relations for 
the soft tri- and bilinear supersymmetry-breaking parameters Aq, Bq such as those inspired 
by specific supergravity scenarios, a framework we term the VCMSSM [7]. Additionally, 
if one assumes universality between mo and the gravitino mass, as in minimal supergravity 
(mSUGRA), the gravitino might be the LSP and constitute the dark matter [8], a framework 
known as the GDM [9-11]. 

It is well known that predicting the masses of supersymmetric particles using precision 
low-energy data is more difficult than it was for the top quark or even the Higgs boson. 
This is because the Standard Model (SM) is renormalizable, so decoupling theorems imply 
that many low-energy observables are insensitive to heavy sparticles [12]. On the other 
hand, supersymmetry may provide an important contribution to loop effects that are rare or 
forbidden within the Standard Model. In fact, we found previously [1] that present data on 
the electroweak precision observables M]y and sin^ 6'eff, as well as the loop induced quantities 
{g — 2)^ and BR(6 s'-f) (sec Ref. [13] for a review), may already be providing interesting 
indirect information on the scale of supersymmetry breaking, at least within the context of 
the CMSSM with a neutralino LSP. In that framework, the range of mo is very restricted by 
the cold dark matter density Q^h"^ determined by WMAP and other observations, for any set 
of assumed values of tan/3,mi/2 and the trilincar soft supersymmetry-breaking parameter 
Aq [14,15]. We fixed mg so as to satisfy this density constraint, 0.094 < r2cDM^^ < 0.129 [16], 
and then analyzed the indirect information as a function of mi/2 for tan/? = 10,50. This 
was done for various discrete values of Aq and as a scan in the (mi/2, ^o) plane. 

Within the CMSSM and using the (then) preferred range nit = 178.0 ±4.3 GcV [17], we 
found previously [1,2] a preference for low values of mi/2, particularly for tan/j = 10, that 
exhibited also a moderate sensitivity to ^40. Our first step in this paper is to update our 
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previous analysis, taking into account the newer preferred range rrit — 172.7 ±2.9 GeV [18], 
and providing a vade mecum for understanding the imphcations of any further evolution in 
the preferred range and experimental error of rrit. The new experimental value of rrit reduces 
substantially the mass expected for the hghtest MSSM Higgs boson, M^, for any given values 
of mi/2, mo, tan /3 and Aq, strengthening the constraints on mi/2- We therefore improve our 
analysis by incorporating the full likelihood information provided by the final results of the 
LEP search for a Standard Model-like Higgs boson [19,20]. 

Other recent analyses [21] in the framework of the CMSSM differ from our analysis by 
the omission of certain observables such as Mw, sin^ ^eff or Mh, or in their treatment of the 
95% C.L. exclusion bound for M^. The other analyses find a preference for somewhat larger 
tan /3, mostly due to the fact that Mw and sin^ ^eff are either ignored or treated differently. 

The main purpose of the present paper is to analyze the sensitivity of the preference for 
a low value of mi/2 to some of the restrictive assumptions we introduced into the analysis, 
exploring the ranges of parameters that would be preferred in alternative NUHM, VCMSSM 
and GDM scenarios. 

The NUHM has two additional parameters as compared to the CMSSM, namely the de- 
grees of non-universality of the soft supersymmetry-breaking scalar masses for the two Higgs 
doublets [6] . They can be traded for two quantities measurable at low energies, such as the 
Higgs mixing parameter and the CT'-odd Higgs boson mass. Ma- We explore here the 
possible sensitivities to these parameters within the NUHM. It would take prohibitive effort 
to analyze systematically all this multi-dimensional parameter space. Therefore, we focus 
here on analyzing a limited number of NUHM scenarios, corresponding to two-dimensional 
subspaces of parameters that generahze specific favoured CMSSM scenarios, with the idea of 
exploring whether the dependences on the additional NUHM variables are capable of mod- 
ifying significantly the CMSSM preference for relatively small values of mi/2 and exploring 
possible preferences for the values of other model parameters. 

On the other hand, in very constrained variants of the MSSM (VCMSSM) in which one 

postulates a relation between the tri- and bilinear soft supersymmetry-breaking parameters: 

Aq — Bo+m.0,^ motivated by simple supergravity, the dimensionality of the model parameter 

space is reduced compared with that in the CMSSM. The supersymmetric vacuum conditions 

then fix the ratio of MSSM Higgs vacuum expectation values tan/5 as a function of m.1/2, mo 

and Aq [7]. We study the cases Ao/mo — 0, 0.75, 3 — VS and 2, which are compatible with 

"'^Our notation for the ^0 and Bq parameters follows that which is standard in supergravity models (see 
e.g. Ref. [3]), namely the coupling in the scalar potential is given by Aq g'^^^ + Bq g^"^^ for the tri- and bi-linear 
superpotential terms g^^^ and g^'^\ respectively. This differs from the sign convention used in many publicly 
available codes, see e.g. Ref. [22]. 
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neutralino dark matter for extended ranges of mi/2, and we discuss the preferred ranges of 
mi/2 and tan (5 in each case. 

In general, yet another relevant parameter, namely the gravitino mass, must be taken 
into account, leading to the possibility that the LSP is the gravitino, in which case it would 
provide dark matter, the GDM scenario. In order to simphfy the analysis of GDM in a 
motivated manner, we restrict our attention to scenarios inspired by minimal supergravity 
(mSUGRA), in which the gravitino mass is constrained to equal mo at the input GUT scale, 
and the trilinear and bilinear soft supersymmetry-breaking parameters are again related by 
Aq — Bq-\- tuq. In the cases we analyze in this paper, namely AQ/niQ = 0, 3/4, 3 — -\/3, 2, the 
regions ^ of the (mi/2,mo) plane allowed by cosmological constraints then take the form of 
wedges located at smaller values of mo than those allowed in CMSSM scenarios [9, 10]. We 
scan here some of the GDM wedges allowed by cosmology, exploring whether the new ranges 
of mo may lead to preferences for different values of mi/2- 

We have performed fits in all scenarios, and our main results are as follows. Within 
the CMSSM, we find that the new, lower value of m^ and new treatment of the constraints 
from the LEP Higgs search do not change greatly the values of mi/2 that were preferred 
previously [1,2]. For example, the 90% C.L. upper bound on mi/2 that we obtain for 
tan /3 = 10 is shifted slightly upwards by about 50 GeV. The minimum value of for 
the global fit is increased, however, primarily because of the increased impact of the LEP 
Mh constraint on the CMSSM parameter space. The tension between Mh and the precision 
electroweak observables would become severe for m* < 170 GeV. The minimum values of 
X^ for tan (5 — and 50 are now very similar. We find that the minimum x^ values remain 
approximately the same also for the intermediate values tan /3 = 20 and tan (3 — 35. On the 
other hand, the upper limit on mi/2 could be increased by as much as about 20% by possible 
future changes in the preferred central value of rrit and likely reductions in its error (assuming 
that the experimental results and theoretical predictions for the precision observables are 
otherwise unchanged), but remains relatively small, in general. 

Within the NUHM, we find that the minimal x^ values are smaller than those for CMSSM 
points with the same value of mi/2, and that x^ is relatively insensitive to Ma but may 
decrease or increase as is varied. The preference for small mi/2 is preserved in at least the 
sparse NUHM sample studied here. However, we do find that mo may differ significantly 
from its preferred range in the CMSSM. Likewise, significantly different values of ji and Ma 
are also possible. In general, within the NUHM scenarios studied, the prospects for observing 

^Thc case Ao = 3 — is motivated by the simplest Polonyi model of Planck-scale supersymmetry 
breaking [23]. 
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sparticles at the LHC or the ILC are similar to those in the CMSSM case, except that in 
some cases the fi may be rather heavier than the x5- 

In most of the VCMSSM scenarios with neutrahno dark matter (NDM), looking along 
the coannihilation strip compatible with WMAP and other cosmological data, we find that 
the preference for small mi/2 noted previously within the CMSSM framework is repeated 
(offering good detection prospects for the LHC and the ILC) , and becomes a preference for 
medium values of tan/3. In addition, there is a tendency for tan/3 to increase with mi/2- 
On the other hand, for Ao/mo = we find larger values of mi/2 at the minimum (which 
is significantly larger than for larger values of Ao/mo), and smaller values of tan/3 which 
are rather constant with respect to mi/2- When Ao/mo — 2, we also observe that there 
are WMAP-compatible VCMSSM models at mi/2 ~ 140 GeV and mo ~ 600 GeV [24] with 
tan /3 ~ 37 that have even lower x^. These occur in the light Higgs funnel, when 2m^o ^ M^, 
and offer some prospects for detection at the Tevatron. 

The preference for small mi/2 and a medium range of tan/3 is also maintained within 
the VCMSSM with the supplementary mSUGRA relation 771,3/2 = mo when the dark matter 
is composed of gravitinos (GDM) and the next-to-lightest supersymmetric particle (NLSP) 
is the fi- In this scenario, the global that is somewhat smaller than along the WMAP 
strips in the VCMSSM with neutrahno dark matter. The prospects for sparticle detection at 
the LHC and ILC are rather similar to those in the previous VCMSSM NDM scenarios, but 
the hght Higgs funnel disappears, reducing the prospects for the Tevatron. We recall that 
the NLSP is metastable in such GDM scenarios, suggesting that novel detection strategies 
should be explored at the LHC and the ILC [25] . 

2 Current experimental data 

In this Section we review briefiy the experimental data set that has been used for the fits. 
We focus on parameter points that yield the correct value of the cold dark matter density, 
0.094 < QcBuh^ < 0.129 [16], which is, however, not included in the fit itself. The data 
set furthermore comprises the following observables: the mass of the W boson, Mw, the 
effective leptonic weak mixing angle, sin^ ^eff , the anomalous magnetic moment of the muon, 
{g — 2)^, the radiative 5-decay branching ratio BR(fo — > 57), and the lightest MSSM Higgs 
boson mass, M^. A detailed description of the first four observables can be found in [1, 13]. 
We hmit ourselves here to recaUing the current precision of the experimental results and 
the theoretical predictions. The experimental values of these obervables have not changed 
significantly compared to [1, 13], and neither have the theoretical calulations. As already 
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commented, due to the new, lower experimental value of m^, it is necessary to include the 
most complete experimental information about Mh into the fit. Accordingly, we give below 
details about the inclusion of Mh and the evaluation of the corresponding values obtained 
from the direct searches for a Standard Model (SM) Higgs boson at LEP [19]. 

In the following, we refer to the theoretical uncertainties from unknown higher-order 
corrections as 'intrinsic' theoretical uncertainties and to the uncertainties induced by the 
experimental errors of the input parameters as 'parametric' theoretical uncertainties. We do 
not discuss here the theoretical uncertainties in the renormalization-group running between 
the high-scale input parameters and the weak scale: see Ref. [26] for a recent discussion in the 
context of calculations of the cold dark matter density. At present, these uncertainties are less 
important than the experimental and theoretical uncertainties in the precision observables. 

Assuming that the five observables listed above are uncorrelated, a fit has been per- 
formed with 



Here denotes the experimental central value of the nth observable {Mw, sin^ 9^^, 

{g — 2)^ and BR(6 — > 57)), R^^"^" is the corresponding CMSSM prediction and (T„ denotes 
the combined error, as specified below. y^M^ denotes the contribution coming from the 
lightest MSSM Higgs boson mass as described below. 

2.1 The W boson mass 

The W boson mass can be evaluated from 



are summarized in the quantity Ar [27]. The prediction for Mw within the Standard Model 
(SM) or the MSSM is obtained by evaluating Ar in these models and solving (j21) in an 
iterative way. 

We include the complete one-loop result in the MSSM [28,29] as well as higher-order QCD 
corrections of SM type that are of 0{aas) [30,31] and [32,33]. Furthermore, we 

incorporate supersymmetric corrections of 0{aas) [34] and of 0{af) [35,36] to the quantity 



A re-evaluation of Mw is currently under way [37] . Preliminary results show good agreement with the 
values used here. 




(1) 




(2) 



where a is the fine structure constant and Gp the Fermi constant. The radiative corrections 
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The remaining intrinsic theoretical uncertainty in the prediction for Mw within the 
MSSM is still significantly larger than in the SM. It has been estimated as [36] 

^^intr,current < 9 MeV , (3) 

depending on the mass scale of the supersymmetric particles. The parametric uncertainties 
are dominated by the experimental error of the top-quark mass and the hadronic contribu- 
tion to the shift in the fine structure constant. Their current errors induce the following 
parametric uncertainties [13,38] 

5m™* = 2.9 GeV aM^,"''"''""™"^"* ^ 17.5 MeV, (4) 

5(A<^S^'=°*) = 36 X 10-^ AMPf''"'^"'^^'^''^™* ^ 6.5 MeV . (5) 

The present experimental value of Mw is [39,40] 

The experimental and theoretical errors for are added in quadrature in our analysis. 

2.2 The effective leptonic weak mixing angle 

The effective leptonic weak mixing angle at the Z boson peak can be written as 
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sin^^efr = 7 fl-Re^) , (7) 

where v^s and agfj denote the effective vector and axial couplings of the Z boson to charged 
leptons. Our theoretical prediction for sin^ ^efr contains the same class of higher-order cor- 
rections as described in Sect. 12.11 

In the MSSM, the remaining intrinsic theoretical uncertainty in the prediction for sin^ 9es 
has been estimated as [36] 

A. 2 /lintr, current ^ n— 5 /d\ 

sm ' < 7 X 10 °, (8) 

depending on the supersymmetry mass scale. The current experimental errors of rrit and 
AcKhad induce the following parametric uncertainties 

(5m™* = 2.9 GeV A sin^ 0^™*'™* ^ 10 x 10-^ (9) 

5(A«™*) = 36 X 10-^ A Sin^ ^P^-.^-had, current _ ^3 ^ ^^q) 

The experimental value is [39,40] 

sin^ ^-P'™"-t = 0.23153 ± 0.00016 . (11) 
The experimental and theoretical errors for sin^ ^etr are added in quadrature in our analysis. 
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2.3 The anomalous magnetic moment of the muon 

The SM prediction for the anomafous magnetic moment of the muon (see [41,42] for reviews) 
depends on the evaluation of QED contributions (see [43] for a recent update), the hadronic 
vacuum polarization and light-by-light (LBL) contributions. The former have been evaluated 
in [44-47] and the latter in [48-51]. The evaluations of the hadronic vacuum polarization 
contributions using e^e~ and r decay data give somewhat different results. In view of the 
additional uncertainties associated with the isospin transformation from r decay, we use here 
the latest estimate based on e~^e~ data [52]: 

af''° = (11 659 182.8 ± 6.3had ± 3.5lbl ± 0.3qed+ew) x 10"!°, (12) 

where the source of each error is labelled. We note that new e^e" data sets have recently 
been published in [53-55], but not yet used in an updated estimate of — Their inclusion 
is not expected to alter substantially the estimate given in (fT^ . 

The result for the SM prediction is to be compared with the final result of the Brookhaven 
{g — 2)^ experiment E821 [56,57], namely: 

a;fP = (11 659 208.0 ± 5.8) x 10"^°, (13) 

leading to an estimated discrepancy 

«)r^ - = (25-2 ± 9.2) X 10"^°, (14) 

equivalent to a 2.7 a effect. While it would be premature to regard this deviation as a 
firm evidence for new physics, it does indicate a preference for a non-zero supersymmetric 
contribution. 

Concerning the MSSM contribution, the complete one-loop result was evaluated a decade 
ago [58]. It indicates that variants of the MSSM with yU < are already very challenged by 
the present data on a^, whether one uses either the e~^e~ or r decay data, so we restict our 
attention in this paper to models with /i > 0. In addition to the full one-loop contributions, 
the leading QED two- loop corrections have also been evaluated [59]. Further corrections at 
the two-loop level have been obtained recently [60,61], leading to corrections to the one-loop 
result that are ~ 10%. These corrections are taken into account in our analysis according to 
the approximate formulae given in [60,61]. 

2.4 The decay b ^ 

Since this decay occurs at the loop level in the SM, the MSSM contribution might a priori be 
of similar magnitude. A recent theoretical estimate of the SM contribution to the branching 
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ratio is [62] 



BR(6 -> S7) = (3.70 ± 0.46) x 10 
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(15) 



where the calculations have been carried out completely to NLO in the MS renormalization 
scheme [63-65], and the error is dominated by higher-order QCD uncertainties. We record, 
however, that the error estimate for BR(6 57) is still under debate, see also Refs. [66,67]. 



For comparison, the present experimental value estimated by the Heavy Flavour Aver- 
aging Group (HFAG) is [68] 



where the error includes an uncertainty due to the decay spectrum, as well as the statistical 
error. The good agreement between ()16|) and the SM calculation (fT^ imposes important 
constraints on the MSSM. 

Our numerical results have been derived with the BR(6 —>■ 57) evaluation provided in 
Ref. [69], which has been checked against other approaches [64,65,70,71]. For the current 
theoretical uncertainty of the MSSM prediction for BR(6 57) we use the value in p5|) . 
We add the theory and experimental errors in quadrature. 

We have not included the decay Bs fi^fi~ in our fit, in the absence of an experimental 
likelihood function and a suitable estimate of the theoretical error. However, it is known 
that the present experimental upper limit: BR(i?s —>■ fi~^fi~) < 2 x 10"'' [72] may become 
important for tan/5 > 40 in the MSSM [73,74]. We mention below some specific instances 
where the decay Bg jX^jj,' may already constrain the parameter space studied [75], and 
note that [1] gives a detailed analysis of its possible future significance. 

2.5 The lightest MSSM Higgs boson mass 

The mass of the lightest CP-even MSSM Higgs boson can be predicted in terms of the other 
CMSSM parameters. At the tree level, the two CP-even Higgs boson masses are obtained as 
functions of M^, the CP-odd Higgs boson mass Ma, and tan [3. For the theoretical prediction 
of Mh we employ the Feynman-diagrammatic method, using the code FeynHiggs [76,77], 
which includes all numerically relevant known higher-order corrections. The status of the 
incorporated results can be summarized as follows. For the one- loop part, the complete 
result within the MSSM is known [78-80]. Computation of the two-loop effects is quite 
advanced: see Ref. [81] and references therein. These include the strong corrections at 
0{atas) and Yukawa corrections at 0{a^) to the dominant one-loop 0{at) term, and the 
strong corrections from the bottom/sbottom sector at 0{ab<ys)- In the case of the b/b sector 



BR(6 57) = {S.Sdtoii) X 10 
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(16) 
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corrections, an all-order resummation of the tan/3 -enhanced terms, tan/5)"), is also 

known [82,83]. Most recently, the 0{atab) and 0{al) corrections have been derived [84] ^. 
The current intrinsic error of Mh due to unknown higher-order corrections has been estimated 
to be [13,81,86,87] 

^^intr.current ^ 3 QeV . (17) 

We show in Fig. C]the predictions for in the CMSSM for tan/3 = 10 (left) and tan/3 = 50 
(right) along the strips allowed by WMAP and other cosmological data [14]. We note that 
the predicted values of depend significantly on Aq. Also shown in Fig. ^ is the present 
95% CL. exclusion limit for a SM-like Higgs boson is 114.4 GeV [19] and a hypothetical 
LHC measurement at Mh = 116.4 ± 0.2 GeV. 

It should be noted that, for the unconstrained MSSM with small values of Ma and values 
of tan /3 that are not too small, a significant suppression of the hZZ coupling can occur in 
the MSSM compared to the SM value, in which case the experimental lower bound on Mh 
may be more than 20 GeV below the SM value [20]. However, we have checked that within 
the CMSSM and the other models studied in this paper, the hZZ coupling is always very 
close to the SM value. Accordingly, the bounds from the SM Higgs search at LEP [19] can 
be taken over directly (see e.g. Refs. [88,89]). It is clear that low values of mi/2, especially 
for tan/3 = 10, are challenged by the LEP exclusion bounds. This is essentially because 
the leading supersymmetric radiative corrections to Mh are proportional to ln(mi/2/'n^t), 
so that a reduction in rrit must be compensated by an increase in mi/2 for the same value 
oiMh. 

In our previous analysis, we simply applied a cut-off on Mh, considering only parameter 
choices for which FeynHiggs gave Mh > 113.0 GeV. However, now that the Mh constraint 
assumes greater importance, here we use more completely the likelihood information avail- 
able from LEP. Accordingly, we evaluate as follows the Mh contribution to the overall 
function ^. Our starting points are the CLs{Mh) values provided by the final LEP results 
on the SM Higgs boson search, see Fig. 9 in [19] ^. We obtain by inversion from CLs{Mh) 
the corresponding value of ^^{Mh) determined from [90] 

ierfc(^^x'(M,)) = CL,(M,) , (18) 

and note the fact that CL,{Mh = 116.4 GeV) = 0.5 implies that x^(116.4 GeV) = as is 

"'A two- loop effective potential calculation has been presented in [85], but no public code based on this 
result is currently available. 

^We thank P. Bechtle and K. Desch for detailed discussions and explanations. 
^We thank A. Read for providing us with the CLs values. 
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Figure 1: The CMSSM predictions for as functions of mi/2 with (a) tan/? = 10 and 
(b) taii/S = 50 for various Aq. A hypothetical LHC measurement is shown, namely = 
116.4 ± 0.2 GeV, as well as the present 95% C.L. exclusion limit of 114.4 GeV. 



appropriate for a one-sided limit. Correspondingly we set x^{Mh > 116.4 GeV) = 0. The 
theory uncertainty is included by convolving the likelihood function associated with x^(M/j) 
and a Gaussian function, normalized to unity and centred around M^, whose width is 

1.5 GeV: 

x\Mh) = -2 log (^J°° e-^'(^)/2 _ _ (19) 

In this way, a theoretical uncertainty of up to 3 GeV is assigned for ~ 95% of all values 
corresponding to one parameter point. The final Xm^ is then obtained as 

x'm, = x'(M,)-x'(116.4GeV) for M, < 116.4 GeV , (20) 

X^, = for M/, > 116.4 GeV , (21) 

and is then combined with the corresponding quantities for the other observables we consider, 
see eq. (H)). 

3 Updated CMSSM analysis 

As already mentioned, in our previous analysis of the CMSSM [1] we used the range = 
178.0 ± 4.3 GeV that was then preferred by direct measurements [17]. The preferred range 
evolved subsequently to 172.7 ± 2.9 GeV [18]. In view of this past evolution and possible 
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future developments, in this Section we first analyze the current situation in some detail, 
emphasizing some new aspects related to the lower value of rrit, and then provide a guide to 
possible future developments. 

The effects of the lower value are threefold. First, it drives the SM prediction of 
Myy and sin^ 9es slightly further away from the current experimental value (whereas {g — 
2)^ and BR(6 57) are little affected). This increases the favoured magnitude of the 
supersymmetric contributions, i.e., it effectively lowers the preferred supersymmetric mass 
scale. Secondly, the predicted value of the lightest Higgs boson mass in the MSSM is lowered 
by the new nit value, see, e.g., Ref. [91] and Fig. H The effects on the electroweak precision 
observables of the downward shift in Mh are minimal, but the LEP Higgs bounds [19, 20] 
now impose a more important constraint on the MSSM parameter space, notably on mi/2- 
In our previous analysis, we rejected all parameter points for which FeynHiggs yielded 
Mh < 113 GeV. The best fit values in Ref. [1] corresponded to relatively small values of Mh, 
a feature that is even more pronounced for the new nit value. Thirdly, the focus-point region 
of the CMSSM parameter space now appears at considerably lower rriQ than previously, 
increasing its importance for the analysis. 

In view of all these effects, we now update our previous analysis of the phenomenological 
constraints on the supersymmetric mass scale mi/2 in the CMSSM using the new, lower 
value of nit and including a contribution from Mh, evaluated as discussed in the previous 
Section. As in Ref. [1] we use the experimental information on the cold dark matter density 
from WMAP and other observations to reduce the dimensionality of the CMSSM parameter 
space. In the parameter region considered in our analysis we find an acceptable dark matter 
relic density along coannihilation strips, in the Higgs funnel region and in the focus-point 
region. We comment below on the behaviours of the function in each of these regions. 

As seen in the first panel of Fig. |2l which displays the behaviour of the function out to 
the tips of typical WMAP coannihilation strips, the qualitative feature observed in Ref. [1] 
of a pronounced minimum in at nii/2 for tan/5 = 10 is also present for the new value of 
nit. However, the curve now depends more strongly on the value of Aq, corresponding to 
its strong impact on Mh, as seen in Fig. [H Values of Ao/mi/2 < —1 are disfavoured at the 
90% C.L., essentially because of their lower Mh values, but Ao/mi/2 = 2 and 1 give equally 
good fits and descriptions of the data. The old best fit point in Ref. [1] had Ao/mi/2 = — 1, 
but there all Ao/mi/2 gave a similarly good description of the experimental data. The 
minimum x^ value is slightly below 3. This is somewhat higher than the result in Ref. [1], 
but still represents a good overall fit to the experimental data. The rise in the minimum 
^See also Ref. [2], where a lower bound of > 111.4 GeV has been used. 
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value of x^; compared to Ref. [1], is essentially a consequence of the lower experimental 
central value of rrit, and the consequent greater impact of the LEP constraint on Mh [19,20]. 
In the cases of the observables Mw and sin^^es, a smaller value of mt induces a preference 
for a smaller value of mi/2, but the opposite is true for the Higgs mass bound. The rise in the 
minimum value of reflects the correspondingly increased tension between the electroweak 
precision observables and the Mh constraint. 

A breakdown of the contributions to from the different observables can be found for 
some example points in Table ^ The best-fit points for tan f3 = 10 and 50 are shown in 
the first and third lines, respectively. The second line shows a point near the tip of the 
WMAP coannihilation strip for tan/5 = 10, the fourth line shows a point at the tip of the 
rapid-annihilation Higgs funnel for tan j3 = 50. The fifth till the seventh row show points 
in the focus point region (see below) for tan f3 = 50 with low, intermediate and high mi/2- 
It is instructive to compare the contributions to at the best-fit points with those at the 
coannihilation, Higgs funnel and focus points. One can see that, for large mi/2 values in 
all the different regions, {g — 2)^ always gives the dominant contribution. However, with 
the new lower experimental value of nit also Mw and sin^ O^s give substantial contributions, 
adding up to more than 50% of the {g — 2)^ contribution at the coannihilation and Higgs 
funnel points. On the other hand, M^ and BR(6 — » 57) make negligible contributions to 

at these points. As seen from the last lines of the Table, the situation may be different 
in the focus-point region for low vnxj^'- the first example given yields a reasonably good 
description of My^, sin^6'efr and even [g — 2)^, while the largest contribution to arises 
from BR(6 57) ^ . This smoothly changes to the behavior for large mi/2 as described 
above also in the focus-point region, as can be seen from the last two rows in Tab. H 

The remaining panels of Fig. El update our previous analyses [1] of the functions for 
various sparticle masses within the CMSSM, namely the lightest neutralino Xi? the second- 
lightest neutralino X2 ^'^d ^^le (almost degenerate) lighter chargino Xa-, the lightest slepton 
which is the lighter stau fi, the lighter stop squark ti, and the gluino g. Reflecting the 
behaviour of the global x^ function in the first panel of Fig. |2l the changes in the optimal 
values of the sparticle masses are not large. The 90% CL. upper bounds on the particle 
masses are nearly unchanged compared to the results for = 178.0 ± 4.3 GeV given in 
Ref. [1]. 

The corresponding results for WMAP strips in the coannihilation, Higgs funnel and focus- 

^We note that, particularly in view of the current uncertainties on mt and m?, and the corresponding 
uncertainties in Ma, the upper limit on the BR(_Bs — > [i^ [i") currently imposes a weaker constraint on the 
CMSSM parameter space than does h 57, even for tan/3 = 50 [74]. 
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320 


90 
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best fit 


2.55 


1.01 


0.12 


0.63 


0.23 


0.52 


10 


880 


270 


1760 


bad fit 


9.71 


2.29 


1.28 


6.14 


0.01 





50 


570 


390 


-570 


best fit 


2.79 


1.44 


0.31 


0.08 


0.91 


0.04 


50 


1910 


1500 


-1910 


bad fit 


9.61 


2.21 


1.11 


6.29 


0.01 





50 


250 


1320 


-250 


focus 


7.34 


0.89 


0.15 


1.69 


3.76 


0.84 


50 


330 


1640 


-330 


focus 


6.06 


1.24 


0.28 


3.21 


1.33 





50 


800 


2970 


-800 


focus 


8.73 


1.92 


0.72 


6.05 


0.04 






Table 1: Breakdown of contributions from the different precision observables to Xtot /^'^ 
some example points. All masses are in GeV. The first and third rows are the best fits for 
tan /5 = 10 and 50, the second row is representative of the coannihilation strip, the fourth 
row is representative of the Higgs funnel region, and the last three rows are representatives 
of the focus point-region. 

point regions for the case tan /3 = 50 are shown in Fig. El The spread of points with identical 
values of Aq at large mi/2 is due to the broadening and bifurcation of the WMAP strip in 
the Higgs funnel region, and the higher set of curves originate in the focus-point region, 
as discussed in more detail below. We see in panel (a) that the minimum value of for 
the fit with rrit = 172.7 ± 2.9 GeV is larger by about a unit than in our previous analysis 
with rrit = 178.0 ± 4.3 GeV. Because of the rise in for the tan/5 = 10 case, however, the 
minimum values of are now very similar for the two values of tan (3 shown here. The dip 
in the function for tan/? = 50 is somewhat steeper than in the previous analysis, since 
the high values of mi/2 are slightly more disfavoured due to their M]^ and sin^ ^gfr values. 
The best fit values of mi/2 are very similar to their previous values. The preferred values 
of the sparticle masses are shown in the remaining panels of Fig. El Due to the somewhat 
steeper x^ behavior, the preferred ranges have slightly lower masses than in Ref. [1]. 

We now return to one novel feature as compared to Ref. [1], namely the appearance of a 
group of points with moderately high x^ that have relatively small mi/2 ~ 200 GeV. These 
points have relatively large values of rriQ, as refiected in the relatively large values of m^^ 
and m^-^ seen in panels (d) and (e) of Fig. El These points are located in the focus-point 
region of the (mi/2,mo) plane [92], where the LSP has a larger Higgsino content, whose 
enhanced annihilation rate brings the relic density down into the range allowed by WMAP. 
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Figure 2: The combined likelihood function for the electroweak observables M-^y, sin^^eff; 
{g~2)^, BR(6 S7), and Mh evaluated m the CMSSM for tan /3 = 10, rrit = 172.7±2.9 GeV 
and various discrete values of Aq, with ttiq then chosen to yield the central value of the relic 
neutralino density indicated by WMAP and other observations. We display the function 
for (a) mi/2, (b) m^o (c) m^o m^±, (d) nif^, (e) nii^ and (f) m-g. 
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Figure 3: As in Fig. \^ but now for tan/5 = 50. 
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By comparison with our previous analysis, the focus-point region appears at considerably 
lower values of mo, because of the reduction in the central value of rrit. This focus-point strip 
extends to larger values of tuq and hence mi/2 that are not shown. The least-disfavoured 
focus points have a Ax^ of at least 3.3 (see the discussion of Table [T] above), and most of 
them are excluded at the 90% C.L. 

Taken at face value, the preferred ranges for the sparticle masses shown in Figs. |21 and El 
are quite encouraging for both the LHC and the ILC. The gluino and squarks lie comfortably 
within the early LHC discovery range, and several electroweakly-interacting sparticles would 
be accessible to ILC(500) (the ILC running at = 500 GeV). The best-fit CMSSM point 
is quite similar to the benchmark point SPSla [93] (which is close to point B of Ref. [94]) 
which has been shown to offer good experimental prospects for both the LHC and ILC [95]. 
The prospects for sparticle detection are also quite good in the least-disfavoured part of the 
focus-point region for tan/3 = 50 shown in Figs. El with the exception of the relatively heavy 
squarks. 

As indicated in Tab. ^above, the minimum values of are 2.5 for tan (3 = 10 and 2.8 for 
tan/5 = 50, found for mi/2 ~ 320,570 GeV and Aq = +mi/2, —mi/2, respectively, revealing 
no preference for either large or small tan/5 ^. We display in Fig. I^the functions for two 
intermediate values of tan/5 = 20, 35, for the cases Aq = 0, ±mi/2- The minima of are 2.2 
and 2.5, respectively, which are not significantly different from the values when tan j3 = 10, 50. 
Thus, this analysis reveals no preference for intermediate values of tan/?, either. The 
minima are found for Aq = 0, — mi/2, respectively. They appear when mi/2 ~ 400, 500 GeV, 
values intermediate between the locations of the minima for tan/5 = 10,50, demonstrating 
the general stability of this analysis. 

In view of the possible future evolution of both the central value of rrit and its experimental 

uncertainty 6mt, we have analyzed the behaviour of the global function for 166 GeV < 

rrit < 179 GeV and 1.5 GeV < 6mt < 3.0 GeV for the case of tan/? = 10 (assuming that the 

experimental results and theoretical predictions for the precision observables are otherwise 

unchanged), as seen in the left panel of Fig. El We see that the minimum value of is 

almost independent of the uncertainty 6mt, but increases noticeably as the assumed central 

value of rrit decreases. This effect is not strong when rrit decreases from 178.0 GeV to 

172.7 GeV, but does become significant for rrit < 170 GeV. This effect is not independent 

of the known preference of the ensemble of precision electroweak data for rut ~ 175 GeV 

within the SM [39,40], to which the observables M^r and sin^ 9^^ used here make important 

^In our previous analysis, we found a slight preference for tan/3 — 10 over tan/3 — 50. This preference 
has now been counterbalanced by the increased pressure exerted by the Higgs mass constraint. 



16 



14 



12 ■ 



10 ■ 



■o 8 
o 



CMSSM, fi>0, m,= 172.";^ 
. tanp = 20, Aj, = 
. tanp = 20, Aj, = +m,,2 
o tanp = 20, A„ = -m,,^ 



' If 



14 



12 ■ 



10 ■ 



■D 8 
o 



CMSSM, fi>0, m,= 172.";^ 
. tanp = 35, A„ = 
. tanp = 35, A„ = +m^^^ 
o tanp = 35, A,, = -m^^ 



••• *, 



200 



400 600 

[GeV] 



800 



1000 



200 



400 600 

[GeV] 



800 



1000 



Figure 4: The functions for tan (3 = 20, 35 and Aq = 0, ±mi/2. 
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90% C.L. upper limit for mi/2 on rrit and its experimental error 6mt, keeping the experimental 
values and theoretical predictions for the other precision observables unchanged. 



contributions. On the other hand, as aheady commented, within the CMSSM there is the 
additional effect that the best fit values of mi/2 for very low rrit result in Mh values that 
are excluded by the LEP Higgs searches [19,20] and have a very large Xm^; resulting in an 
increase of the lowest possible value for a given top-quark mass value. This effect also 
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Table 2: Breakdown of contributions from the different precision observables to x^^^ for 
some example points with rrit = 168,173,178 GeV, Srrit = 2.5 GcV and tan/3 = 10. All 
masses shown are in GeV. The fifth column shows the value for the corresponding point, 
and the last column shows the contribution of this Mh value. The values of Aq were 
selected so as to minimize Xtot /^^^ corresponding value of mi/2- 



18 



increases the value of mi/2 where the function is minimized. This is analyzed in more 
detail in Table El where we show the breakdown of the different contributions to ioT 
rrit = 168, 173, 178 GeV for 6mt = 2.5 GeV and tan/? = 10. The Aq values are chosen so 
as to minimize Xtot each choice of mi/2. For = 168 GeV, Xtot exhibits only a shallow 
and relatively high minimum, and and BR(6 — + s'j) give the largest contribution for 
low mi/2, shifting smoothly to large contributions from Mw, sin^ 9es and {g — 2)^ for larger 
mi/2- For nit = 173 GeV, a more pronounced minimum of Xtot appears for relatively small 
mi/2 values. For lower mi/2, again and BR(6 — > 57) give large contributions, whereas for 
higher values this shifts again to M^r, sin^ 9^^ and — after passing through a minimum 
with a very good fit quality where no single contribution exceeds unity. The same trend, just 
slightly more pronounced, can be observed for rrit = 178 GeV. Finally, in the right panel of 
Fig. El we demonstrate that the 90% C.L. upper limit on mi/2 shows only a small variation, 
less than 10% for rrit in the preferred range above 170 GeV Finally we note that the 
upper limit on mi/2 is essentially independent of 6mt 

It is striking that the preference noted earlier for relatively low values of mi/2 remains al- 
most unaltered after the change in rrit and the change in the treatment of the LEP lower limit 
on Mh- There seems to be little chance at present of evading the preference for small mi/2 
hinted by the present measurements of Mw, sin^ 9cs, BR(6 — > S7) and ((7 — 2)^, at least within 
the CMSSM framework. It should be noted that the preference for a relatively low SUSY 
scale is correlated with the top mass value lying in the interval 170 GeV ^ < 180 GeV. 

4 NUHM Analysis 

In the NUHM, one may parametrize the soft supersymmetry-breaking contributions to the 
squared masses of the two Higgs multiplets, follows: 

= ml{l + 6i) , (22) 

where ml is the (supposedly) universal soft supersymmetry-breaking squared mass for the 

squarks and sleptons. As already mentioned, the increase of the dimensionality of the NUHM 

parameter space compared to the CMSSM, due to the appearance of the two new parameters 

61^2, makes a systematic survey quite involved. Here, as illustrations of what may happen 

in the NUHM, we analyze some specific parameter planes that generalize certain specific 

^"^The plot has been obtained by putting a smooth polynomial through the otherwise slightly irregular 
points. 

^^Note added: this analysis demonstrates, in particular, that incorporating the latest global fit value 
mt = 172.5 ± 2.3 GeV [96] would have a negligible effect on our analysis. 
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CMSSM points. We note that certain combinations of input parameter choices lead to soft 
SUSY-breaking Higgs mass squares which are negative at the GUT scale. When either 
+ /i^ < or ml + /i^ < 0, the point is excluded, so as to ensure vacuum stability at the 
GUT scale [6]. 

Since it is the value of mi/2 that affects most importantly the masses of the sparticles that 
might be observable at the LHC or ILC, our primary objective is to investigate whether the 
introduction of extra NUHM parameters affects significantly the preference for small mi/2 
found previously within the CMSSM: see Figs. El and [7| After satisfying ourselves on this 
point, we subsequently explore the possible dependences on and /i: see Fig. |H1 In order 
to present our results we use parameter planes with generic points that do not necessarily 
satisfy the CDM constraint. Exhibiting full parameter planes rather than just the regions 
where the neutralino relic density respects the WMAP limits (we indicate these strips in 
the plots) provides a better understanding of the dependences of the function on the 
different NUHM variables. It also provides a context for understanding the branchings of 
the function visible in Fig. El which are due to the bifurcations of the WMAP strips in 
the parameter planes. We also note that, in NUHM models with a light gravitino where the 
CDM constraint does not apply, generic regions of these parameter planes may be consistent 
with cosmology. 

In view of our primary objective. Fig. IHlshows two examples of (mi/2, ""^o) planes for fixed 
values of /i. Ma, Aq and tan (3 (top row) and two examples of (^«i/2, ^) planes for fixed values 
of mo. Ma, Aq and tan (3 (bottom row). In both the two top panels, the left boundaries of the 
shaded regions are provided by the LEP lower limit on the chargino mass, the upper bounds 
on mi/2 are provided by the GUT stability constraint, and the lower edges of the shaded 
regions are provided by the stau LSP constraint. The colour codings are as follows. In each 
panel, the best fit NUHM point that respects the WMAP constraints on the relic neutralino 
density is marked by a (blue) plus sign, and the (blue) cross indicates the CMSSM values of 
(mi/2,'Tio) [or {iTLi/2,f^)] for the chosen values of the other parameters. The green (medium 
grey) regions have Ax^ < 1 relative to the minimum when the WMAP/CDM constraint is 
not employed. Hence, some points in this region may have a lower than our best fit point 
when the CDM constraint is employed. 

In all four panels of Fig. IHl our best CDM fit, denoted by the plus sign, is within 1 sigma 
of the overall minimum x^, and hence lies within the green region. The yellow (light grey) 
regions have A^^ < 3.84, and the black points have larger values of Ax^ relative to the 
absolute minimum. Traversing the regions with Ax^ < 1, 3.84, there are thin, darker shaded 
strips where the relic neutralino density lies within the range favoured by WMAP. That is. 
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in these regions, is within 1 or 3.84 of the minimum when the WMAP/CDM bound 
is included. The blue cross must always lie within these regions. Our sampling procedure 
causes these WMAP strips to appear intermittent. In the top right panel of Fig. IHl we 
note two vertical tramlines, which are due to rapid annihilation via the direct-channel A 
pole. Since Ma is fixed in each of these panels, there is always a value of mi/2 such that 
2m^o ^ Ma, in principle even for tan/5 = 10. We note that the analogous tramlines are 
invisible in panel (a), because they have a Ax^ > 4 and thus would be located in the black 
shaded region. 

In the lower two panels, large values of fi are excluded due to the GUT constraint, large 
values of mi/2 are excluded by the stau LSP constraint, and low values of fi and mi/2 are 
exuded by the chargino mass limit. In the lower left panel, large values of mi/2 are excluded 
because the fi becomes the LSP, whereas in the right panel our computation was limited 
to mi/2 < 1000 GeV, thus producing the right-hand boundary. Within the regions allowed 
by these constraints, the same colour codings are used. In the lower right panel, one sees 
clearly the effect of the pseudoscalar funnel at mi/2 ~ 680 GeV. In the lower left panel, this 
possibility is excluded by the GUT stability constraint. 

The planes in Fig. IHl have been defined such that the CMSSM points marked by (blue) 
crosses in the different panels of Fig. IHllie at the minima of the CMSSM functions shown 
in Figs. El and El They enable us to study whether the CMSSM preference for relatively 
small mi/2 may be perturbed by generalizing to the larger NUHM parameter space. In each 
case, we see that the CMSSM point lies close to the best NUHM fit, whose is lower by 
just 0.00, 0.02, 0.72 and 0.08, respectively. We also note that the ranges of mi/2 favoured at 
this level are quite close to the CMSSM values. Thus, in these cases, the introduction of two 
extra parameters in the NUHM does not modify the preference for relatively small values of 
mi/2 observed previously in the CMSSM. In the top left panel for tan/5 = 10, we see that 
the preferred range of mo is also very close to the CMSSM value. On the other hand, we 
see in the top right panel that rather larger values of mo would be allowed for tan /5 = 50 at 
the Ax"^ < 1 level. This is due to the insensitivity of the annihilation cross section to mo in 
the funnel due to rapid annihilation via the pseudoscalar Higgs boson A. We also see in the 
bottom two panels that quite wide ranges of fi would be allowed for either value of tan/5 

Fig. El displays four analogous NUHM planes, specified this time by values of /x, and 

Aq in the top row and mo. Ma, Aq in the bottom row that do not correspond to minima of 

the function for the CMSSM with the corresponding values of tan/5. These examples 

-"^^In all panels of Fig. IHl the assumed values of Ma are sufficiently large that Bg /^^M^ currently does 
not impose any useful constraint [75] . 
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Figure 6: Sample NUHM scenarios shown in the (1711/2,1^0) plane (top row) and 
(mi/2,fJ') plane (bottom row). The CMSSM points shown in the left (right) column cor- 
respond to the best fit points for tan/3 = 10 (50). The other parameters are given in the 
plots. The green [medium grey] (yellow [light grey]) regions have Ax^ < 1(3.84), whilst the 
black regions have larger Ax^ . The strips where the neutralino relic density respects the 
WMAP limits have darker shadings. The blue plus sign marks the best-fit NUHM point that 
respects the relic density bounds, and the CMSSM point is marked with a blue cross. 
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Figure 7: Additional sample NUHM scenarios shown in the (mi/2,mo) plane (top row) and 
(mi/2,fJ') plane (bottom row). The colour coding is the same as in Fig.{^ 
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were studied in detail in [6] , and enable us to explore whether there may be good NUHM fits 
that are not closely related to the best CMSSM fits. In the top panels, the left boundaries 
are due to the chargino constraint, and the bottom boundaries are due to the stau LSP 
constraint. In the left panel, the right boundary is due to GUT stability, but in the right 
panel it is due to a sampling limit. In the bottom left panel the GUT, stau and chargino 
constraints operate similarly as in Fig. IHl and the tail at low fi and large mi/2 is truncated 
by the GUT stability constraint. In the bottom right panel, the top boundary is due to 
GUT stability, the bottom boundary to the stau, and the boundary at large mi/2 is another 
sampling limitation Within the allowed regions of Fig. the colour codings are the same 
as in Fig. IHl The best fit CDM point lies within the Ax^ < 1 green regions in the top left 
and bottom right panels, whereas in the upper right panel the best fit point has Ax^ slightly 
larger than 1, and its Ax^ is even greater in the bottom left panel. 

In the {mi/2, mo) planes shown in the top row, we see that the ranges of mi j 2 favoured at 
the Ax^ < 1 level are again limited to values close to the best-fit CDM values. The range of 
mo for a given Ax"^ is somewhat restricted for tan/3 = 10 (top left), but is again considerably 
larger for tan/9 = 50 (top right). As for the (mi/2,yu) planes in the bottom row, we see in 
the left panel for tan/3 = 10 that the range ofmi/2 is again restricted at the Ax^ < 1 level, 
whereas the range of /i is almost completely unrestricted. A similar conclusion holds in the 
bottom right panel for tan/3 = 50, though here the range of mi/2 is somewhat broader 

Having established that the CMSSM preference for small values of mi/2 is generally 

preserved in the NUHM, whereas different values of mo and are not necessarily disfavoured, 

we now study further the sensitivity to /i and M4 via the four examples of (/i. Ma) planes 

shown in Fig. |H1 In each case, we have made specific choices of Aq, tan/3, mi/2 and mo. In 

the two panels on the left, these correspond to the best CMSSM fit along the corresponding 

WMAP strip. The examples on the right were studied in [6]. In each case, we restrict our 

attention to the regions of the plane that have no vacuum instability below the GUT scale. 

This constraint provides the near-vertical right-hand edges of the coloured regions, whereas 

the other boundaries are due to various phenomenological constraints. The near-vertical 

boundaries at small in the top panels are due to the LEP chargino exclusion, and those 

in the bottom panels are due to the stau LSP constraint. The boundary at low M^ in the 

top left panel is also due to the stau LSP constraint, whereas that in the top right panel is 

again the GUT stability constraint. 

^•^We note that, in this example, the CMSSM point is excluded by the stau LSP constraint. 
^"^In all panels of Fig. [7| the assumed values of Ma are again sufhciently large that Bg fJ-^fJ-^ currently 
does not impose any useful constraint [75] . 
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Within the allowed regions of Fig. |H1 the colour codings are the same as in Fig. IHl We see 
that in the top left panel the WMAP strip runs parallel to the lower boundary defined by 
the stau LSP constraint. The best fit NUHM point has = 1.19, which is somewhat less 
than two units smaller than for the CMSSM point. This is hardly significant, and suggests 
that the absolute minimum of the NUHM lies at a similar value of mi/2- As seen from 
the location and shape of the green region with Ax^ < 1, the fit is relatively insensitive to 
the magnitudes of /i and M4, as long as they are roughly proportional, but small values 
of ^/Ma are disfavoured. In contrast, for the larger value of mi/2 shown in the top right 
panel of Fig. |H1 we see that low values of /i are preferred. However, the minimum value of 

= 5.12 in the NUHM is not much lower than in the CMSSM, even though it occurs for 
significantly smaller values of both yU and M4 

Turning now to the bottom left panel of Fig. |H1 for tan /? = 50 and Aq = 580 GeV, with 
mi/2 and ttiq again chosen so as to minimize (i-e., to reproduce the corresponding best-fit 
point), we note several features familiar from the two previous panels. The WMAP strip 
clings close to the left boundary of the allowed region, apart from an intermittent funnel 
straddling the M4 = 2m-j^o line. The minimum of = 1-13 for the NUHM is somewhat 
smaller than for the CMSSM. The Ax^ < 1 NUHM region is a lobe extending away from the 
origin at small fi and M^. Similar features are seen in the bottom right panel for tan/5 = 50 
and Aq = 0, except that the Ax^ < 1 lobe extends up to rather larger values of fi and Ma 

These examples show that, although the absolute values of /i and Ma are typically 
relatively unconstrained in the NUHM their values tend to be correlated, often with 
a restricted range for their ratio: M^/yU ~ 1.4, < 1,~ 2 at the Ax^ < 1 level in the first 
three panels of Fig. |H1 On the other hand, the correlation in the fourth panel takes the form 
Ma ~ i(/i-400 GeV). 

To conclude this Section, we make some remarks about the preferred masses of sparticles 
and their possible detectability within the NUHM framework, in the light of the above x^ 
analysis. Since the ranges of mi/2 favoured within the CMSSM are also favoured in the 
NUHM, one should expect that the LHC prospects for detecting the gluino and several other 
sparticles may also be quite good in the NUHM. On the other hand, the greater uncertainties 

""^^We recall that, in this case, the NUHM WMAP strip has two near-horizontal branches straddling the 
Ma = 2mp^o contour, with the upper branch heading to large Ma at small /i, features not seen clearly in 
this panel because of the coarse parameter sampling. 

^^We note, however, that the lower ranges of Ma ^ 300 GeV in the two bottom panels of Fig.|Slare likely 
to be excluded by the current upper limit on BR(_Bs fi^ fj,~) [75], once the experimental likelihood is made 
available and combined with the corresponding theoretical errors. 

^^The prospects for an indirect determination of Ma and /i using future Higgs-sector measurements have 
been discussed in [97]. 
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in mo,/U and Ma suggest that the prospects for sparticle studies at the ILC may be more 
variable within the NUHM. These remarks are borne out by Figs. El and EB which display 
functions for various sparticle masses in a selection of NUHM scenarios. Fig. 121 presents 
masses in the four NUHM scenarios shown in Fig. |H1 in which the CMSSM points correspond 
to the best-fit points from Sect. El and Fig. ^1 presents masses in two of the scenarios shown 
in Fig. IHl 

In each panel of Fig. El we display the functions for the masses of the Xi, ')<^/Xii ^i, 
and for NUHM parameters along the WMAP strips in the corresponding panels of Fig. El 
Since there are several branches of the WMAP strips in some cases, the functions are 
sometimes multivalued. In the top left panel of Fig. El we see well-defined preferred values 
for the sparticle masses, with the gluino and stop masses falling comfortably within reach 
of the LHC, and the Xi^^i and possibly also the X2 ^'^^ xf within reach of the ILC(500). 
When ~ 4.5, new branches of the x^ function appear, corresponding to a branching of 
the WMAP strip around a rapid-annihilation funnel when Ma = 559 GeV ~ 2m^o. This 
funnel is not visible in Fig. El but would appear in the black-spotted region of large A^^. 
The ILC(IOOO) would have a good chance to see even the lighter stop. Turning to the top 
right panel of Fig. El we see that the branching due to the rapid-annihilation funnel appears 
at much lower Ax^, reflecting the closeness of the funnel to the best-fit point in the top right 
panel of Fig. El In this case, whereas the g should be observable at the LHC, the ti might 
well be problematic The Xi would be kinematically accessible at the ILC(500), but the 
fi might well be too heavy: the rises in the branches of its x^ function at larger masses 
reflect the extension of the WMAP strip to large uiq that is seen in the corresponding panel 
of Fig. El In this particular scenario, the X2 xf would probably not be observable at the 
ILC(500). The ILC(IOOO) on the other hand, would have a high potential to detect them. 
The bottom left panel of Fig. El has the most canonical x^ functions: the gluino and stop 
would very probably lie within reach of the LHC and the fi within reach of the ILC(500), 
whereas the X2 ^"^^ Xi might be more problematic. Again the ILC(IOOO) offers much better 
opportunities here, possibly even for the lighter stop. Finally, the prospective observabilities 
in the bottom right scenario would be rather similar to those in the top right scenario: we 
again see that, as one moves away from the coannihilation strip, the fi may become much 
heavier than the x\i ^ind too heavy to observe at the ILC(500). The ILC(IOOO) should, on 
the other hand, offers very good prospects. 

Fig. ^1 presents a similar analysis of sparticle masses of the two favoured scenarios in 
Fig. El namely in the two left-hand panels. In these cases, we show the variations of the x^ 
""^^We recall that it is thought to be observable at the LHC if it weighs less than about 1 TeV. 



27 




10r 



T3 

o 



- f •; 



NUHM, X _,. = 1.83 



2 - 



mm w 



m„ = 100, = 559, A,, = 600, tanp = 10 

_l , I I I I I , I I , I I , I L 

500 1000 1500 

mass [GeV] 



10 



CO 
T3 
O 



2 - 



NUHM, = 2.71 



./ . /. 



/ 



2000 



m„ = 390, = 583, A,, = -580, tan|3 = 50 

_l , I I \ I I , I I , I I , I L 

500 1000 1500 

mass [GeV] 



2000 



Figure 9: The likelihood along the WMAP strips in the sample NUHM scenarios shown 
in Fig. O as a function of the masses of the Xi, xVxf) ^i; ^.''^d g. The branchings in the 
curves reflect the corresponding branchings in the WMAP strips in Fig. O 



28 



functions for the different masses as one follows the WMAP strip to larger values of Ma- In 
the left panel of Fig. we display the masses of the X2 ^"^^ Xi (which are nearly equal) 
in black, the mass of the Xs in pink, the masses of the X4 s'^nd xt (which are nearly equal) 
in red, the mass of the ti in yellow (with black border), and in blue. In each case, 
the + sign of the same colour represents the best fit in the CMSSM for the same values of 
'^i/2;'^0;^o and tan/5. The fact that the minima of the NUHM lie somewhat below the 
CMSSM points reflect the fact that the NUHM offers a slightly better fit, but the difference 
is not significant. In this case, the preferred masses of the X2 and xf are almost identical 
to the best-fit CMSSM values, and the same would be true for the Xi and ^, which are not 
shown. The masses of the Xa? ^1 and A are also very similar to their CMSSM values, but the 
X4 and X2 may be significantly heavier. In addition to the above sparticle masses, the right 
panel also includes the mass of the fi in orange. In this case, whereas the masses of the x? 
(not shown), xH/xi and g (not shown) preferred in the NUHM are similar to their values 
at the best-fit CMSSM point, this is not true for the other sparticles shown. The A boson 
may be considerably lighter, the X35 x\ and the X2' niay be either lighter or heavier, and the 
fi and ti might be significantly heavier for points along the Higgs funnel visible in Fig. |H1 
Thus, in this case the prospects for detecting some sparticles at the LHC or ILC may differ 
substantially in the NUHM from the CMSSM. 

To summarize: these examples demonstrate that, although the preferred value of the 
overall sparticle mass scale set by mi/2 may be quite similar in the NUHM to its CMSSM 
value, the masses of some sparticles in the NUHM may differ significantly from the corre- 
sponding CMSSM values. 

5 VCMSSM Analysis 

As an alternative to the above NUHM generalization of the CMSSM, we now examine par- 
ticular CMSSM models with the additional constraint Bq = Aq — itlq motivated by minimal 
supergravity models, namely the VCMSSM framework introduced earlier. We still assume 
that the gravitino is too heavy to be the LSP. The extra constraint reduces the dimen- 
sionality of the VCMSSM parameter space, as compared with the CMSSM, facilitating its 
exploration. In the CMSSM case, the electroweak vacuum conditions can be used to fix 
and Ma as functions of mi/2,fnQ and Aq for a large range of fixed values of tan/5. On the 
other hand, in the VCMSSM case the expression for Bq in terms of Aq and mo effectively 
yields a relation between |yu| and Ma that is satisfied typically for only one value of tan/5, 
for any fixed set of mi/2,TnQ and Aq values [7,98]. 
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Figure 10: The likelihood along the WMAP strips in the sample NUHM scenarios shown 
in the left panels of Fig. \^ as a function of the masses of the x^/xf (black), Xs (pink), 
xVxf (red), fi (orange) [omitted from the left panel], fi (yellow with black border) and A 
boson (blue). The branchings in the curves in the right panel reflect the corresponding 
branchings in the WMAP strips in the bottom left panel of Fig. The crosses indicate the 
corresponding best fit points in the CMSSM. 



As already mentioned, motivated by {g — 2)^ and (to a lesser extent) BR(6 57), we 
restrict our attention here to the case fi > 0. As is well known, other phenomenological 
constraints tend to favour tan/3 > 5, see e.g. Refs. [91,99]. This condition is generally 
obeyed along the WMAP coannihilation strip for neutralino dark matter in the VCMSSM if 
one assumes Aq > 0, in which case the resultant value of tan (3 tends to increase with mi/2 
and mo along the WMAP strip. We have studied the choices Ao/niQ = 0, 0.75, 3 — ^/3 and 2. 
In this Section we restrict our attention to these cases, and in the next Section we compare 
the VCMSSM results with the corresponding gravitino dark matter scenarios. 

Since in the CMSSM the value of x^ tends first to decrease and then to increase with 
mi/2, but does not vary strongly with tan/?, we would expect the x^ function to exhibit 
a similar dependence on mi/2 also in the VCMSSM scenario. This effect is indeed seen in 
the first panel of Fig. ^2 there are well-defined local minima at mi/2 ~ 400 to 600 GeV, 
as Ao/mo varies from to 2. However, for the latter value of Ao/rrio, we notice some 
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isolated (red) points with mi/2 ~ 140 GeV and much lower ~ 2}^ At these points, 
which barely survive the LEP chargino limit, rapid annihilation through a direct-channel 
light-Higgs pole brings the neutralino relic density down into the WMAP range [24]. The 
remaining panels of Fig. [TT] display the functions for the masses of the Xi; ^i; X2; X? ; 
and g. Their qualitative features are similar to those shown earlier for the CMSSM, with 
the addition of the exceptional low-mass rapid-annihilation points. In these VCMSSM NDM 
scenarios, the LHC has good prospects for the g and ii and the ILC(500) has good prospects 
for the Xi ^^^^ ^i? whereas the prospects for the X2 xf would be dimmer, except at the 
isolated rapid-annihilation points. The ILC(IOOO), on the other hand, would have a good 
chance to detect the X2 the xf? depending somewhat on Ao/mQ. These points might 
also be accessible to the Tevatron, in particular via searches for gluinos. 

We find no analogous focus-point regions in the VCMSSM. When Ao/mQ is large, the 
RGE evolution of fi does not reduce it, even when mo is very large For smaller Ao/rriQ, 
the value of tan/3 fixed by the electroweak vacuum conditions in the VCMSSM becomes 
small: tan /3 < 5 when rriQ is large. In this case, as in the CMSSM, the focus-point region is 
not reached. 

In order to understand better the variation of with mi/2 in Fig. and in particular to 
understand its relatively low value at the rapid light-Higgs annihilation points with mi/2 ~ 
140 GeV [24], we display separately in Fig. [T^lthe dependences of (a) M]v, (b) sin^6'cfT, (c) 
BR(6 57), (d) {g-2)^ and (e) Mh on mi/2 for the case Ao/mo = 2 ^i. Along the VCMSSM 
WMAP strip, we see that Mw prefers a very low value of mi/2, with the rapid-annihilation 
points slightly disfavoured, whereas sin^ 9es prefers a range of somewhat larger values of 
mi/2, with the rapid-annihilation points slightly favoured. However, we then see that both 
BR(6 57) and {g — 2)^ independently strongly disfavour mi/2 ~ 200 GeV, whereas the 
rapid-annihilation points fit these measurements very well. The same tendency is observed 
for Mh. 

These behaviours can be understood by referring to panel (f ) of Fig. El where the regions 

disfavoured by 6 57 and favoured by {g — 2)^ are shaded green and pink (darker and 

lighter grey), respectively. The shaded {g — 2)^ region represents a 2-a deviatation based 

on ()14|) . while the dashed lines represent the region favoured at the 1-a level. The LEP 

Higgs constraint is a diagonal (red) dot-dashed line, while the near- vertical black dashed line 

^^Similar points appear in the CMSSM, but at values of ylo/mi/2 niuch larger than those considered in [1] 
and here. 

20This is true also in the CMSSM. 

^^The values of tan/3 in the VCMSSM are too small for Bg M^A*^ currently to make any significant 
contribution to the fimction [75]. 
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Figure 12: Results as functions ofmi/2 along the VCMSSM WMAP strip for neutralino dark 
matter when Ao/niQ = 2, for (a) Mw, (b) sin^ 9cs, (c) BR(6 — > 57), (d) — and (e) Mh- 
Note the NDM points at low mi/2 ~ 140 GeV that simultaneously fit very well BR(6 S7) 
and {g — 2)^. In panel (f) we display phenomenological constraints in the {mi/2, rrio) plane for 
the VCMSSM with Ao/mo = 2, including both the NDM WMAP strip (blue) and the CDM 
wedge (yellow). The regions disfavoured by b ^ s'-f and favoured by {g — 2)^ are shaded green 
and pink (darker and lighter grey), respectively, the LEP Higgs constraint is a near-vertical 
(red) dashed line, and the (blue) dotted lines are contours o/tan/?, as fixed by the VCMSSM 
vacuum conditions. 
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shows the LEP constraint on the chargino mass. The pale (blue) shaded strip is favoured 
by WMAP for NDM. Below this strip, there is a red shaded region in which the LSP is the 
fi and therefore excluded. Below the fi LSP region, the gravitino is the LSP [9]. In the 
unshaded portion of the GDM region, the next-to-lightest supersymmetric particle (NLSP) 
will decay into a gravitino with unacceptable effects on the abundances of the light elements 
and is excluded by BBN [9,10,100,101]. The pale (yellow) shaded wedge is favoured for 
gravitino dark matter as this region is allowed by BBN constraints. Finally, the black dotted 
curves labeled 20, 25, 30 and 35 correspond to the values of tan (3 required by the VCMSSM 
vacuum conditions. We see that the rapid-annihilation tail of the WMAP strip rises at low 
mi/2 into a region allowed by 6 ^ 57, favoured by {g — 2)^ and tolerated by Mh- It is 
the synchronized non-monotonic behaviour of these last three observables that explains the 
similar non-monotonic behaviour of along the NDM WMAP strip in Fig. ^2 and the low 
value of for the isolated rapid-annihilation point at mi/2 ~ 140 GeV [24]. This is in fact 
the best overall fit point in this VCMSSM scenario, as seen in Fig. ^2 

The preferred ranges of mi/2 seen in Fig. ^2 correspond, through the VCMSSM vacuum 
conditions, to preferred ranges in tan (3. As seen in Fig. El these increase with the chosen 
value of Ao/mo, as does the correlation with mi/2. For Ao/mQ = (top left panel), tan /3 ~ 7, 
increasing to tan/? ~ 10, 15,32 for Aq/uiq = 0.75,3 — -\/3, 3, respectively. In the last case, 
descending the VCMSSM WMAP strip to lower mi/2, whereas we see that exceeds 10 
for mi/2 < 350 GeV, we see again the isolated dark (red) rapid- annihilation points with 
mi/2 ~ 140 GeV [24], which have relatively large tan (3 ~ 37. 

We conclude that the extra constraint imposed in the VCMSSM modifies but does not 
remove the preference found within the CMSSM for small mi/2. Within the VCMSSM 
with neutralino dark matter, the minimum of usually occurs along the generic WMAP 
coannihilation strip at mi/2 ~ 500 GeV. However, when Ao/itlq = 2, we find lower values 
of in the rapid light-Higgs annihilation region with mi/2 ~ 140 GeV. The preferred 
value of tan (3 varies between ~ 7 and ~ 32 on the generic WMAP strip, depending on the 
value of Ao/mo, but tan/? ~ 37 in the light Higgs-pole annihilation region for Ao/mo = 2. 
These points offer prospects for a gluino discovery at the Tevatron: all the other preferred 
parameter sets offer good prospects for observing sparticles at the LHC and ILC(500). 

6 GDM Analysis 

The relation Aq = Bq + mo is just one of the further conditions on supersymmetry-breaking 
parameters that would be imposed in minimal supergravity (mSUGRA) models. The other 
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Figure 13: Illustration of the preferred regions in the space of VCMSSM models for (a) 
Aq = 0, (h) Ao/niQ = 0.75, (c) A^/mQ = 3 — and (d) Ao/mo = 2. In each case, the 
red points show the minimum, the green points have Ax^ < 1, the orange points have 
Ax^ < 3.84, and the black points have larger x^ . 
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is the equality between mo and the gravitino mass. So far, we have imphcitly assumed 
that the gravitino is sufficiently heavy that the LSP is always the lightest neutralino Xi 
and the cosmological constraints on gravitino decays are unimportant. However, this is not 
always the case in mSUGRA models. Indeed, in generic mSUGRA scenarios, as seen in the 
bottom right panel of Fig. ^[ in addition to a WMAP strip where the Xi is the LSP as we 
have assumed so far, there is a wedge of parameter space at lower values of mo (for given 
choices of mi/2 and the other parameters), where the gravitino is the LSP. In this case, there 
are important astrophysical and cosmological constraints on the decays of the long-lived 
NLSP [10, 100, 101], which is generally the lighter stau fi in such mSUGRA scenarios 

Fig. El displays the function for a sampling of GDM scenarios obtained by applying 
the supplementary gravitino mass condition to VCMSSM models for Ao/rrio = 0, 0.75, 3-^3 
and 2, and scanning the GDM wedges at low mo. These wedges are scanned via a series of 
points at fixed (small) mo and increasing mi/2- We note that there is a marginal tendency for 

to increase with increasing mo, though this is not as marked as the tendency to increase 
with mi/2, and that the scan lines are more widely separated for the smaller values of Aq. 
Comparing Figs. ^2 and El we see that lower x^ values may be attained in the GDM cases. 
The third panel of Fig. El and last panel of Fig. El illustrate how this comes about in the 
case Aq/itiq = 2: there is a large contribution to from 6 — > 57 in the NDM for small 
mi/2 that is absent in the GDM, which strongly prefers the combination of smaller mo and 
smaller tan/3 found in the GDM models 

As seen in Fig. the global minimum of for all the VCMSSM GDM models with 
Aq/itiq = 0, 0.75, 3 — and 2 is at mi/2 ~ 450 GeV. However, this minimum is not attained 
for GDM models with larger mo, as they do not reach the low-mi/2 tip of the GDM wedge 
seen, for example, in the last panel of Fig. El In general, we see in the different panels of 
Fig. El that, as in the CMSSM, there are good prospects for observing the g and perhaps 
the ii at the LHC, and that the ILC(500) has good prospects for the Xi and fi, though these 
diminish for larger rriQ. The ILC(IOOO), again, offers much better chances also for large mo. 
We recall that, in these GDM scenarios, the fi is the NLSP, and that the Xi is heavier. The 
fi decays into the gravitino and a r, and is metastable with a lifetime that may be measured 
in hours, days or weeks. Specialized detection strategies for the LHC were discussed in [25]: 
this scenario would offer exciting possibilities near the fi pair-production threshold at the 
ILC. 

^^There are also non-mSUGRA scenarios in which the NLSP is the Xi- Such models are subject to similar 
astrophysical and cosmological constraints, but we do not consider them here. 

^■^The values of tan (3 in these GDM models are also too small for Bg 1^^ ^-^ currently to make any 
significant contribution to the function [75] . 
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Figure 14: The dependence of the function on mi/2 for GDM scenarios with Ao/niQ 
0,0.75,3 - and 2, scanning the regions where the lighter stau Ti is the NLSP, shown 
a function of (a) mi/2, (b) rn^o^, (c) m^o andm^±, (d) rrif-^, (e) nii^, and (f) rrig. 
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As discussed above, a feature of the class of GDM scenarios discussed here is that the 
required value of tan/5 increases with mi/2- Therefore, the preference for relatively small 
mi/2 discussed above maps into an analogous preference for moderate tan/3, as shown in 
Fig. 1151 The different panels are for the four choices Aq/itiq = 0,0.75,3 — -\/3 and 2. In 
each case, the red point indicates the minimum of the function, the green points have 
Ax^ < 1 corresponding to the 68 % confidence level, the orange points have Ax^ < 3.84 
corresponding to the 95 % confidence level, and the black points have larger A^^. We see 
that, at the 95 % confidence level 

300 GeV < mi/2 < 800 GeV, 15 < tan /5 < 27 (23) 

in this mSUGRA class of GDM models. 

7 Conclusions 

Precision electroweak data and rare processes have some sensitivity to the loop corrections 
that might be induced by supersymmetric particles. As we discussed previously in the con- 
text of the CMSSM [1,2], present data exhibit some preference for a relatively low scale of 
soft supersymmetry breaking: mi/2 ~ 300 . . . 600 GeV. This preference is largely driven by 
{g — 2)f^, with some support from measurements of Mw and sin^ 9es- In this paper we have re- 
evaluated this preference, in the light of new measurements of rrit and M\y, and treating more 
completely the information provided by the bound from the LEP direct searches for the Higgs 
boson. The preference for mi/2 ~ 300 . . . 600 GeV is maintained in the CMSSM, and also 
in other scenarios that implement different assumptions for soft supersymmetry breaking. 
These include the less constrained NUHM models in which the soft supersymmetry-breaking 
scalar masses for the two Higgs multiplets are treated as free parameters as well as more con- 
strained VCMSSM models in which the soft trilinear and bilinear supersymmetry-breaking 
parameters are related. The same preference is also maintained in GDM models motivated 
by mSUGRA, where the LSP is the gravitino instead of being a neutralino as assumed in 
the other scenarios. 

Whilst mi/2 is quite constrained in our analysis, there are NUHM scenarios in which itlq 
could be considerably larger than the corresponding values in the CMSSM, and significant 
variations in fi and Ma are also possible. Within the CMSSM and NUHM, we find no 
preference for any particular range of tan f3, but the preferred values of mi/2 in the VCMSSM 
and GDM scenarios studied here correspond to intermediate values of tan/5 ~ 15 to 30. 

The ranges of mi/2 that are preferred would correspond to gluinos and other sparticles 
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Figure 15: Illustration of the preferred regions in the space of mSUGRA-motivated GDM 
models for (a) Aq = 0, (h) Ao/mo = 0.75, (c) Ao/mQ = 3 — and (d) Ao/mo = 2. In each 
case, the red points show the minimum, the green points have Ax^ < 1, the orange points 
have Ax^ < 3.84, and the black points have larger x^- 
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being light enough to be produced readily at the LHC Many sparticles would also be ob- 
servable at the ILC in the preferred CMSSM, VCMSSM and GDM scenarios considered, but 
the larger values of mo allowed in some of the NUHM scenarios would reduce the number 
of sparticle species detectable at the ILC, at least when operated at 500 GeV, whereas the 
ILC at ^/s = 1000 GeV covers the full range for some sparticle species. There are also 
prospects for detecting supersymmetry at the Tevatron in some special VCMSSM models 
with neutralino dark matter. 

We re-emphasize that our analysis depends in considerable part on the estimate of the 
Standard Model contribution to {g — 2)^ based on e+e~ annihilation data, that we assume 
in this paper. Our conclusions would be weakened if the Standard Model calculation were 
to be based on r decay data. Additional e'^e" data are now coming available, and it will 
be important to take into account whatever update of the Standard Model contribution 
to {g — 2)^ they may provide. However, the measurement of Mw is increasing in impor- 
tance, particularly in the light of the recent evolution of the preferred value of rrit. Future 
measurements of and rrit at the Tevatron will be particularly important in this regard. 
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